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Abstract 

The decay of the Isoscalar Giant Monopole Resonance (ISGMR) in nuclei is studied by means 
of a nonlinear classical model consisting of several noninteracting nucleons (particles) moving in 
a potential well with an oscillating nuclear surface (wall). The motion of the nuclear surface is 
described by means of a collective variable which appears explicitly in the Hamiltonian as an 
additional degree of freedom. The total energy of the system is therefore conserved. Although the 
particles do not directly interact with each other, their motions are indirectly coupled by means 
of their interaction with the moving nuclear surface. We consider as free parameters in this model 
the degree of collectivity and the fraction of nucleons that participate to the decay of the collective 
excitation. Specifically, we have calculated the decay width of the ISGMR in the spherical nuclei 
208 Pb, 144 Sm, 116 Sn and 90 Zr. Despite its simplicity and its purely classical nature, the model 
reproduces the trend of the experimental data which show that with increasing mass number the 
decay width decreases. Moreover the experimental results (with the exception of 90 Zr) can be well 
fitted using appropriate values for the free parameters mentioned above. It is also found that these 
values allow for a good description of the experimentally measured 112 Sn and 124 Sn decay widths. 
In addition, we give a prediction for the decay width of the exotic isotope 132 Sn for which there is 
experimental interest. The agreement of our results with the corresponding experimental data for 
medium-heavy nuclei is dictated by the underlying classical mechanics i.e. the behaviour of the 
maximum Lyapunov exponent as a function of the system size. 

PACS numbers: 24.30.Cz,24.60.Lz,05.45.Pq,21.10.Re 
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I. INTRODUCTION 



The study of the decay of nuclear giant resonances (GRs) is a subject of extensive the- 
oretical and experimental investigation GRs having a large degree of collectivity 
are offered for the study of dissipation of collective motion by single particle motion in a 
many body quantum system. Moreover, it is known that regular and chaotic dynamics 
usually coexist in nuclear excitations and chaotic dynamics is expected to dominate in GRs 
% [5| • The latter is due to the time-dependence induced in the nucleons dynamics and to 
the correlations among the nucleons originating from the residual interaction. By studying 
the decay of GRs we can get information about the role that chaotic dynamics plays in 
dissipation of the collective motion. In the damping of the collective motion in GRs several 
mechanisms are at work at the mean field level (one-body dissipation), where the decay is 
into single particle motion (p-h excitation) including particle emission (escape width) and 
beyond (n-body dissipation - spreading width), where the decay occurs through coupling to 
progressively more complicated states starting from 2p-2h excitations (collisional damping). 
The relative contribution of these damping processes is under investigation. In this work we 
consider explicitly one-body dissipation and implicitly part of the n-body dissipation and 
investigate the influence of the underlying classical dynamics on it. 

The effects of the chaotic dynamics on the dissipation of nuclear collective motion and 
the consequent excitation of the nucleons has been considered in a number of previous 
works. We mention for one body dissipation the application of the Fermi accelerator model 
6|, |7[ , the wall formula 8| and its generalizations and computer simulations of classical and 
quantal systems of noninteracting nucleons colliding with a moving boundary undergoing 
periodic adiabatic oscillations 9|. In these earlier works the coupling of the slowly moving 
collective degree of freedom to the fast moving independent particle degrees of freedom has 
been considered disregarding self-consistency and energy conservation. Subsequently, the 
damping of the collective motion coupled self-consistently, ensuring energy conservation, to 
the single-particle motion has been studied using the Vlasov equation 101 ] - To clarify the 
relationship between chaos at the microscopical level and damping of collective motion, a 
classical version of the vibrating potential model ll| for finite nuclei has been considered and 



the case of monopole oscillations along with others has been investigated 



12 



13| . Several 



other calculations of the damping have been carried out later by exploring the phase space 
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[l^ . by considering two-body interactions Q, [3] etc. 

We focus on the decay of the Isoscalar Giant Monopole Resonances (ISGMR) (breathing 
mode, where neutrons and protons oscillate in phase) which are also useful for the extraction 
of compression moduli of nuclei Ka and from them the compressibility of nuclear matter 
\v\ . Recent and ongoing experimental activity at different facilities including the cyclotron 
at Texas A&M and at RCNP at Osaka supplement and extend the previous data clarifying 
the relation between the symmetry energy and Ka- We will consider initially the data for 
the ISGMR of the spherical medium-heavynuclei: 208 Pb, 144 Sm, 116 Sn and 90 Zr Q, 3, EJ 
and those for the 112 Sn and 124 Sn isotopes 2l|. There has been a large number of theoretical 



22, \23 
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From 



approaches for the study of the energy and the decay width of ISGMR 
these studies we mention _self-consistent Hartree-Fock plus Random Phase Approximation 
(RI 



31 



A) with Gogny 25J, |26( or Skyrme 



23 



27, m 



29( or unitary correlation operators 



30] effective interactions and self-consistent relativistic mean field plus RPA calculations 
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33[ as well as calculations that go beyond mean fie 



d approximation, i.e. Extended 



351 ] . In the latter reference, by 



RPA 34j, relativistic RPA with particle-phonon coupling 
extending the covariant density functional theory with the consideration of quasiparticles 
dressed with a cloud of low-lying collective excitations, the influence of the phonon coupling 
terms to the damping of GRs has been studied and led to noticeable fragmentation and 
spreading width. The decay widths extracted from all the above calculations are smaller 
(although to a different extent) from the experimental ones. 

In this work we investigate the role that chaoticity plays in the decay of ISGMR of 
the above mentioned spherical nuclei focusing mainly on one-body processes and using a 
classical model that exhibits chaotic behaviour and that is based on the model of Refs. 



12 



13]. Chaoticity in relation with the decay of the ISGMR has also been considered in 



Refs. 



31 



361 ] in the context of the relativistic mean-field model. 



The considered model is a simple classical Hamiltonian system which consists of particles 
moving in a Woods-Saxon well with an oscillating surface. This oscillation represents the 
ISGMR state in which the nuclear surface is oscillating radially, i.e. the spherical symmetry 
of the nucleus is preserved. The particles can interact with the surface and can exchange 
energy with it. The collective variable which describes the motion of the nuclear surface 
appears explicitly in the Hamiltonian of the system as an additional degree of freedom, 
therefore the total energy of the system is conserved. Moreover, in our study we take into 
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account the escape of nucleons from the oscillating well. Our model therefore takes explicitly 
into account one-body dissipation (Landau damping and particle decay) and partly many- 
body dissipation. Chaoticity is introduced by the time-dependence of the potential and by 
the indirect coupling of the motion of the nucleons by means of their interaction with the 
oscillating nuclear surface, which gives rise to some of the correlations among nucleons that 
are due to the residual interaction. In Refs. 



37 



our model has been used to elucidate 



the relationship of relaxation dynamics with the structure of the phase space of the system. 
We calculate the decay width of the ISGMR in the spherical nuclei 208 Pb, 144 Sm, 116 Sn and 
90 Zr. It is found that under the appropriate assumptions the trend of the experimental 
data can be reproduced, i.e. the decay width is a decreasing function of the mass number A. 
Moreover, with the proper choice of the free parameters, a good quantitative agreement with 
the experimental results can be obtained. This is attributed to the dependence of the mean 
value of the maximum Lyapunov exponent on the size of the nucleus. With the same choice 
of the parameters the decay width of three other isotopes of Sn ( 112 Sn, 124 Sn and 132 Sn) is 
calculated. Finally, the fraction of the width due to the escape is deduced for the above 
nuclei. In Sec. II we present a brief description of the model and the method of calculation. 
In Sec. Ill our results are presented and discussed whereas in Sec. IV the conclusions are 
drawn amd some prospects for further study within our model are mentioned. 

II. DESCRIPTION OF THE MODEL AND METHOD OF CALCULATION 

The model we consider consists of a system of N noninteracting nucleons (particles) of 
mass m moving in a Woods-Saxon well. The nucleons exchange energy with the nuclear 
surface (wall of the potential well), which is considered to move in a harmonic oscillator 



potential with angular frequency a> [12| , Il3l . 1371 . |38| . A mass M is assigned to the nuclear 



surface and its motion around the equilibrium position Rq is described by the collective one- 
dimensional degree of freedom R. Since the potential well is of finite depth, nucleons can 
escape from it yielding a contribution to the decay width. The Hamiltonian of the system 
equals 



N 
i=l 



Pri Vq 



+_ exp (o^) 
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The angular momentum Lr of the collective degree of freedom R is taken equal to zero, i.e. 
the nuclear surface is not rotating. The angular momenta Lj of the particles are constants 
of the motion since the oscillating potential retains its spherical symmetry. The values of 
the parameters we use equal a = 0.67 fm, V = AhMeV and R Q = X.IA^^fm. The above 
values of a and Vq correspond to the nucleus 208 Pb. There is a slight variation of a and Vq 
with the mass number A but for the present study it is neglected. For the angular frequency 
uj we use the frequency of the ISGMR obtained from the experimentally measured centroid 
energy E x . The values of the frequency we use for the nuclei 208 Pb, 144 Sm, 116 Sn and 90 Zr 
equal 13.96Me1/ • k~\ 15A0MeV ■ fT 1 , 15.85MeV ■ h' 1 and ITMMeV ■ h' 1 respectively 
(see Table H| [rsl. 19]). The parameters M/m and N will be considered as free parameters 
and their optimal values will be determined by comparing our results to the experimental 
data. More precisely, free parameters are the quantities M/mA ■ 100% and N/A ■ 100% (see 
EqsJT2land [T3]) . The equations of motion are 



Va 



cxp 



Pr 



Pri 

m 

n-R 

a 



exp 
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mrf 
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y N 



exp (^) 



1 S [l + exp (adBj 



and have been solved using a leapfrog algorithm 



Mu 2 (R-Rq) 



(2) 
(3) 
(4) 
(5) 



391 ] . With a simple rescaling of the equations 
of motion we find that the relevant parameters are N, M/m, R$/a and Vq/ \ma 2 uj 2 ). 

For each nucleus we consider a microcanonical ensemble of 5000 initial conditions selected 
according to the following prescription: 

1. The initial energy Er associated with the collective variable R 

E R 



(6) 



is taken equal to the experimentally measured centroid energy E x of the ISGMR. The 
corresponding initial momentum pr is taken equal to zero, i.e. the harmonic oscillator 
describing the collective motion is initially at one of its extremal points: 



R = Rq 



' 2E X 



2 ' 



(7) 
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2. The initial position rj of the nucleon is uniformly distributed in the interval [0, R]. 
The nucleon is a proton with a probability P = Z/A and a neutron with a probability 
I — P. 



3. The initial relative energy AEi of the nucleon with respect to the bottom of the 
potential well is selected using the density of energy of a Fermi gas, namely 



n{Ei) 



(8) 



2EI 



where Ep is the Fermi energy, which for protons is given by 



E - k 
Ef ^ ~ 2^ 



and for neutrons by 



E 



F(n) 



2m 



3ir 2 



,A 



V 



(9) 



(10) 



The volume V of the system is taken equal to that of a sphere of radius Rq. The initial 
kinetic energy 7$ of the nucleon is therefore given by 

Vo 



Ti = -Vo + AEi + 



exp 



n-R 



(11) 



and its partition in radial and angular part is determined by the angle 0, which is 
uniformly distributed in the interval [0, ir] (0^ is the angle between the vectors of the 
initial position and initial velocity), obtaining the radial component of the momen- 
tum p ri and the angular momentum Lj from the equations p ri = \j2mTi cos 6^ and 



Li = \J ImTiTi sin 9i respectively. We should note that we have also tried other initial 
momentum distributions (i.e. uniform or Boltzmann) and it has been found that the 
results are not sensitive to the particular choice of the initial momentum distribution. 

In our study we neglect the Coulomb barrier. However, we have carried out calculations 
with the Coulomb barrier for the nucleus 208 Pb and our results for the decay width did not 
change appreciably. For each nucleus, we consider as free parameter the quantity 

M 



A 



mA 



■ 100%, 



(12) 



which is the percentage of the mass assigned to the collective degree of freedom and conse- 
quently can be thought of as a measure of the degree of collectivity of the motion. Apart 



from A, we also consider as a free parameter the quantity 

2 = ^ ■ 100%, (13) 

which measures the percentage of nucleons that participate into the damping of the collective 
motion. For each nucleus and for a grid of values of 5 and A, we calculate the mean value 
< R(t) > of the collective variable R over the selected ensemble of N to t = 5000 orbits. We 
also calculate the mean value of the maximum Lyapunov exponent < Ai >, which is the 
most commonly used measure of the intensity of chaotic dynamics, using the same initial 
conditions and integration time as the ones we used in the calculation of < R(t) >. < X\ > 
quantifies how fast initially closed trajectories deviate in phase space and it is related to 
the decay width T (or the relaxation time). The comparison of < Ai > and T reveals the 
microscopic origins of relaxation dynamics as well as the relative importance of the regular 
regions in phase space 33, |38j . In order to estimate the decay width, we take the Fourier 
transform of < R(t) > and calculate the standard deviation a of the corresponding power 
spectrum in a frequency region around uj where the power spectrum has a maximum. In 
order to obtain the width T, which corresponds to the full width at half maximum (FWHM) 



of an equivalent Gaussian fit to the power spectrum, we multiply o by 2\/m4 ~ 2.355. 

It seems that the use of a classical model is quite satisfactory and the violations of Fermi 
statistics are not so important. In particular, the effect of the spin on the momentum and 
energy distribution is not expected to be significant. Results are not sensitive to the type 
of the initial momentum distribution. In fact, the small influence of the Fermi statistics 
can be drawn by the small deviation of our results from the experimental data (see section 
III). Blocki and coworkers 9] have come to the same conclusion by performing computer 
simulations and have shown that results were similar for the excilation of a classical or 
quantal gas of particles by a time-dependent potential well. The model used in the present 
work differs from the one of Baldo et al 12j, [l3| mainly in the consideration of diabatic 
motion for the nuclear surface and realistic geometry for the nuclear potential. Moreover, 
in the present model the number of particles that collide with the wall as well as the mass 
attributed to the oscillating surface have been treated as free parameters. 



S 





FIG. 1: Decay width T of the ISGMR as a function of the parameters S and A for the nuclei (a) 
208 Pb, (b) 144 Sm ,(c) 116 Sn and (d) 90 Zr. In (a), (b) and (c) the lines on which our results coincide 
with the experimental ones are also shown (solid line). The line on which the value of the decay 
width is maximum is shown as a dashed line. 



III. RESULTS AND DISCUSSION 

For the nuclei 208 Pb, 144 Sm, 116 Sn and 90 Zr, and for a 9 x 9 grid of values for S and A, we 
calculate the decay width T and the mean value of the largest Lyapunov exponent < Ai >. 
Our results for T and < Ai > are shown in Figs. [1] and [3] respectively. In Figs. [1] and [2l 
the lines on the H-A plane on which our results coincide with the experimental ones are also 
shown for the case of the nuclei 208 Pb, 144 Sm and 116 Sn. 

We first observe that the relationship between < Ai > and T is in most cases what we 
would expect for a completely chaotic system. More specifically, for all values of S and A, 



9 



1 00 -■ 
90- 




i 1 1 1 1 1 1 1 1 — 

20 40 60 80 100 

E(%) 



FIG. 2: The lines on the S-A plane on which our results for the decay width coincide with the 
experimental ones for the nuclei 208 Pb (solid line, slope 0.63), 144 Sm (dashed line, slope 0.75) and 
116 Sn (dotted lines, slopes 0.51 and 0.46). 

both T and < Ai > decrease as A increases. Moreover, for each of the four nuclei studied, 
as H and A vary, < Ai > and T vary with the same monotonicity except from a region with 



H > 50% and A < 30%. Different monotonicities of < Ai > and T have been reported in 37] 
and can be attributed to the existence of regular regions in the phase space. In all cases, the 
decay width T seems to have its maximum value on a line with slope ~ 1/2 on the 5 — A 
plane (see FigfTJ dashed lines). 

Comparing our results for the decay width with the corresponding experimental values, 
which are shown in Table [J we observe that the trend of the experimental results is repro- 
duced, i.e. the amplitude T decreases with increasing A for all values of H and A. At a 
quantitative level, in the case of the nuclei 208 Pb and 144 Sm, our results coincide with the 
experimental ones on a curve which is close to a straight line (see FigsJHa),(b) and Figj2j). 
The corresponding slopes determined by a least squares fit equal 0.63 and 0.75 respectively. 
Above (below) these lines, our results underestimate (oveserstimate) the experimental data. 
In the case of the nucleus 116 Sn, our results coincide with the corresponding experimental 
ones on two lines with slopes 0.51 and 0.46 (see FigJT](c) and Figf2]). For (5, A) points be- 
tween these lines our results overestimate the experimental ones, whereas for points outside 
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FIG. 3: The mean value < Ai > of the maximum Lyapunov exponent corresponding to the 
calculations of Figfflfor the nuclei (a) 208 Pb, (b) 144 Sm ,(c) 116 Sn and (d) 90 Zr. 



this region the experimental data are underestimated. In the case of the nucleus 90 Zr, our 
results for T are smaller than the corresponding experimental ones for all values of S and A. 
It seems that our model is more successful in the description of the decay width of medium- 
heavy nuclei. This may be due to the decreased collectivity of the ISGMR in lighter nuclei 
(A < 100) as well as to the increasing importance of higher-body processes in these nuclei. 
In the considered model, as it has already been mentioned, although the nucleons do not 
interact with each other, their motions are indirectly coupled by means of their interaction 
with the moving nuclear surface and therefore part and not the whole effect of the higher 
body processes is taken into account. 

We must mention that we have also considered the case of 1 particle moving inside the 
potential well, as in Ref. 37|. In that case, although the trend of the experimental data is 
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Nucleus 


E x (MeV) 


r exp (MeV) 




208p b 


13.96 ± 0.20 [181 


2.88 ±0.20 [18] 


0.11 + 0.042 [401 


144 Sm 


15.40 ± 0.30 [181 


3.40 ±0.20 [18] 




116 Sn 


15.85 ± 0.20 [181 


5.27 ±0.25 [18] 




90 Zr 


17.81 + 0.20, -0.32 [12] 


7.86 + 0.89,-1.41 [19] 


~ 0.08 [41] 


112 Sn 


15.67 ±0.11 [21] 


5.18 + 0.40,-0.04 [21] 




124 Sn 


15.34 ±0.13 [21] 


5.00 + 0.53,-0.03 [21] 





TABLE I: Experimental results for the excitation energy E x , the decay width T exp and the fraction 



^ esc of the width due to nucleon escape for the ISGMR in the nuclei 2U °P 



112 Sn and 124 Sn (Results for 116 Sn coincide with those reported recently 



144 



Sm. 



116 



Sn, 9U Zr, 



reproduced (r decreases with increasing A), the results for T are very small compared to 
the experimental ones. Therefore, the consideration of many particles turned out necessary 
for the quantitative description of the decay widths. 

In order to determine a region in the (5, A) plane where there is an optimal agreement 
between our results and the corresponding experimental values, we assume that all four 
nuclei can be described using the same values of A and H, i.e. all the nuclei considered 
exhibit the same degree of collectivity and the same fraction of nucleons that take part in 
the damping. We calculate the sum £ 2 of the squares of the relative differences between our 
results and the experimental data 



£ 2 



Tpj, — rp b ( cxp ) \ / r Sm — r 5m ( exp ) \ / — r5 n ( exp ) \ / r Zr — r^ r ( Cxp ) 
rpb( exp ) j \ r,s m ( CX p) j \ rs- n ( cxp ) J \ r^ r ( Cxp ) 



(14) 

as a function of the parameters S and A (Sn = 116 Sn). The results are shown in Fig. [4] 
From this figure we observe that we have a good agreement with the experimental results in 
a region on the (H, A) plane which is around a line with a slope approximately equal to 0.69. 
The optimal agreement with the experimental results is obtained where £ 2 has its minimum 
value, i.e. for H ~ 69% and A ~ 62%. We should note that by assuming a uniform sphere 
oscillating with linear radial displacement and velocity field, an oscillating mass M equal to 
3/5 (0.6) of the total mass of the nucleus is deduced (scaling model) 44], |45(. For the above 
mentioned values of H and A, the calculated decay widths for the nuclei 208 Pb, 144 Sm, 116 Sn 
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and 90 Zr equal 3.08MeV, 3.94MeV, 4.64MeV and 5.05MeV respectively. The above results 
along with the corresponding experimental ones are shown in Figj5](a). We observe that the 
discrepancy between our results and the experimental data increases as A decreases. This 
is possibly due to our assumption that all nuclei exhibit the same degree of collectivity, 
which should be modified when considering lighter nuclei. We should also mention the role 
of the n-body processes and the gradual change of the shape of the ISGMR response which 
is observed for A < 100 losing its symmetric single peaked Gaussian-like character. As 
for the dependence of T on A in our model, it is dictated by the classical mechanics of 
the system since the mean maximum Lyapunov exponent < Ai > of the system has the 
same dependence on A, as can be seen from FigjS^b). Comparing our results for the width 
T of ISGMR with those from other approaches mentioned in the Introduction, we realize 
that they are generally smaller (although to a different extent) than ours. For example, 
recent results obtained with self-consistent Hartree-Fock plus RPA approach with SkM* 



effective interaction and fully taking into account the continuum [461 ] yield the following 
results for the width: r( 208 Pb) = 2.5MeV and r( 90 Zr) = 3.0MeV [47j, whereas the values 
of the width T obtained with the relativistic quasiparticle RPA along the lines of Ref. 48] 
using the effective interaction DDME2 derived from an effective Lagrangian with density- 
dependent meson-nucleon vertex function are r( 208 Pb) = 1.54MeV and r( 90 Zr) = 2.18MeV 



49j . We realize that both sets of values of the width are in qualitative agreement with the 



experimental ones but they are smaller from the latter and from the values derived with 
our classical model. This is expected since only lp-lh space has been considered in RPA. 
Recently, as mentioned in the introduction, the influence of phonon coupling terms to the 
damping of Giant Resonances has been studied and led to a noticeable fragmentation and 



increase of the width 



35j (for instance for 208 Pb, T equals 2 MeV and 3 MeV with Relativistic 



RPA and Relativistic RPA plus Particle- Vibration Coupling respectively). 

For the above values of H and A (H = 69% and A = 62%), we have also calculated the 
decay widths of the ISGMR of two other Sn isotopes, namely 112 Sn and 124 Sn. Our calcula- 
tion yields for the widths 4.68 MeV and 4.05 MeV respectively, whereas the experimentally 



obtained values equal 5.18 MeV and 5.00 MeV respectively [21] (Recently new results for 
the centroid energy E x and width have been reported 43j . The discrepancy with the values 
of E x and V for 112 Sn in Tabled has been restored; see also conclusions). In addition, for 
the same values of H and A, we made a prediction for the decay width of the exotic nu- 
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FIG. 4: The sum £ 2 of the squares of the relative differences between our results and the experi- 
mental data for the nuclei 208 Pb, 144 Sm, 116 Sn and 90 Zr (see EqEJ). 



cleus 132 Sn for which there is increased experimental interest. Our calculation yields T equal 
to 3.87 MeV. In our calculation, we have used as centroid energy E x the value 15.29MeV, 
which is calculated theoretically in Ref. 26( by performing Hartree-Fock plus Random Phase 
Approximation calculation using the D1S parametrization of the Gogny two-body effective 
interaction. We should mention that results for the decay width Y of the isotopes 112 Sn, 124 Sn 
and 132 Sn have been derived from other theoretical approaches (see i.e. 26|, |35|, |48j, l50l. |51|). 
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means of the two approaches reported 



491 ] . Recently, a decay width of 132 Sn 



The values for the decay width of 132 Sn derived 
above equal 3.21MeV and 2.13MeV respectively 
equal to 3.09MeV has been reported in 35J by fitting the theoretical strength distribution 
with a Lorentzian. These results are smaller than the value found in our calculations. 

We have also calculated the fraction /i esc of the total decay width which is due to the 
escape of nucleons. It is defined as 

r-r' 



(15) 



where T is the total decay width and T' is the decay width without taking into account the 
escape of particles. The latter is extracted from the mean value < R(t) > of the collective 
variable over the ensemble of orbits that do not lead to nucleon escape, i.e. (R(t)) = 
J2 Rj{t) /N to t where Rj(t) is the time series of the collective variable in the j—th orbit and 

3 



14 



10- 
9- 
8-\ 
7- 
6- 
5- 
4- 
3H 



(a) 



5 



80 100 120 140 160 180 200 220 

A 



> 

CD 



4.0 

3.8H 

3.6 

3.4 

3.2 



2.8 



2.6 



(b) 



80 100 120 140 160 180 200 220 

A 



FIG. 5: (a) The calculated values of the total decay width T for the nuclei 208 Pb, 144 Sm, 116 Sn and 
90 Zr (circles) for the values of S and A (S — 69%, A ~ 62%) that minimize the quantity £f (se e 
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Eq ll4p . The corresponding experimental results (squares) with their errors are also shown 
(b) The corresponding values of the mean maximum Lyapunov exponent < Ai >. 



the summation is performed only over the orbits in which all the nucleons remain confined 
in the potential well at time t. In the orbits where a nucleon escapes from the potential well, 
the motion of the escaping nucleon is followed until its distance from the nucleus reaches a 
cutoff value, which has been taken much larger than the radius of the nucleus. Then, the 
two degrees of freedom describing the motion of the escaping nucleon are considered to be 
decoupled from the motion of the nucleons confined in the potential well, i.e. the escaping 
nucleon moves away from the nucleus with constant energy and it is therefore excluded 
from the sum of eq 5. Once the nucleon escapes and travels at constant energy away from 
the nucleus, it carries energy away from the system, thus contributing to the dissipation 
of the collective degree of freedom. The results are shown in Figj6l From this Figure, it 
can be seen that in the region of S and A values in which agreement of our results for T 
with the experimental data is optimal, \x esc is between 0.20 and 0.28 for all nuclei under 
consideration. These values are larger than the correspon ding exp erimental ones which are 
available only for the nuclei 208 Pb and 90 Zr (see Table [J 
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531 ] . Inclusion of 

a Coulomb barrier leads to a decrease of /i esc by a fraction ranging from « 21% (for 90 Zr) 
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FIG. 6: The fraction fi esc of the total decay width due to the escape of nucleons as a function of 
the parameters H and A for the nuclei (a) 208 Pb, (b) 144 Sm ,(c) 116 Sn (d) 90 Zr. 

to m 35% (for 208 Pb). However, as our calculation is purely classical it does not account 
for the quantum tunneling of particles through the potential barrier. One should add that 
the experimental values for \i esc reported above refer to the direct escape width whereas our 
model probably includes part of the statistical decay width. The partial widths for direct 
neutron escape of the ISGMR of 90 Zr, 124 Sn and 208 Pb have also been calculated recently 



4l| within a continuum extended RPA approach with the use of a phenomenological mean 



field, the Landau-Migdal p-h interaction and some partial self-consistent conditions. 
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IV. CONCLUSIONS AND PROSPECTS 



In conclusion, in this work we studied the decay of the Isoscalar Giant Monopole Reso- 
nance (ISGMR) for several spherical nuclei using a classical model, in which a number of 
noninteracting nucleons (particles) move in an oscillating potential well and can exchange 
energy with the nuclear surface (wall of the potential well). The motion of the nuclear 
surface is described by a collective variable which appears explicitly in the Hamiltonian of 
the system as an additional degree of freedom. The total energy of the system is therefore 
conserved. Since there is no interaction between the particles, the model takes explicitly 
into account only one-body processes. Nevertheless, the particles interact indirectly with 
each other by means of the coupling of their motion to the motion of the oscillating well and 
therefore dissipation beyond the one-body is expected to be taken into account. Moreover, 
for medium-heavy nuclei, with the proper choice of the free parameters of the system (the 
degree of collectivity and the fraction of nucleons that participate to the decay), our results 
show a good agreement with the experimental data, at least comparable to the agreement of 
more sophisticated quantum models. The dependence of the decay width of the ISGMR on 
the size of the nucleus has been found to be related to the characteristics of the underlying 
classical dynamics, i.e. to the maximum Lyapunov exponent. Its value is determined by the 
presence of chaotic regions in the phase space. Chaoticity is included in the dynamics of 
GRs due to the time-dependence and to the residual interaction. We expect that at least 
qualitatively, our results should hold in the quantum systems. Recently, there have been 
experimental results for the ISGMR of several isotopes of Sn from RCNP in Osaka 431 ] . The 
theoretical calculations of the energy give values that are larger than the experimental ones. 
These imply larger values of the compression moduli. Work is in progress to investigate 
with our model in more detail the underlying classical dynamics of ISGMR in these isotopes 
that differ in the symmetry energy and predict their decay widths. We should add that our 
model properly adjusted can be used for the study of the decay of ISGMR of other nuclei 
including those far from stability as well as for the study of the Isovector Giant Monopole 
Resonances. 
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